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AN APPROXIMATE SOLUTION FOR THE RADAR ECHO PULSE RESPONSE 
FOR PLANETARY RADAR ALTIMETERS 

By W. Thomas Bundick 
Langley Research Center 

SUMMARY 

An approximate solution for the radar echo pulse response has been developed to 
compute the expected return pulse amplitude and shape for planetary radar altimeters. 
Results obtained with the solution show excellent agreement with results obtained by 
numerical integration; in addition, this present solution requires much less computer 
time. The new solution is applicable over a wider range of altitudes and pulse widths than 
the solution previously presented by Harrington and Stanley in NASA TN D-5220; the pre- 
sent solution is also more accurate at low altitudes and for long pulse widths. 

INTRODUCTION 

For pulse-type radar altimeters, the arrival time of the beginning of the received- 
pulse leading edge is a measure of the shortest distance to the surface. This time of 
arrival, however, is difficult to measure because most pulse- radar-altimeter range 
trackers measure the arrival time of some other point on the received pulse, frequently 
some point on the leading edge. However, the shape and rise time of the leading edge is a 
function of the surface backscatter characteristics, which vary with incidence angle and 
terrain. These variations in the shape of the return pulse produce errors in the altitude 
measurements known as terrain bias errors. 

Knowledge of the radar echo pulse response as a function of altitude and surface 
backscatter characteristics can be used to determine variations in the altimeter terrain 
bias error. In reference 1 the radar pulse response was derived using the model of ref- 
erence 2 for the surface backscatter characteristics. The solution, as shown in equa- 
tion (C4) of reference 1 and in equation (A34) of appendix A of this report, is an approxi- 
mation which is valid for certain ranges of normalized time and altitude. In some cases 
it is desirable to know the response under conditions which violate the assumptions made 
in reference 1. Therefore, the purpose of this report is to offer a more general solution 
which is less restrictive than that of reference 1, as well as to make a comparison 
between predicted radar return waveforms for the two solutions. 



For the convenience of the reader, a slightly modified version of Harrington and 
Stanley's derivation of the radar echo step response (ref. 1) is presented in appendix A of 
this report, and a review of Muhleman's model for the radar backscatter function (ref. 2) 
is presented in appendix B. 


SYMBOLS 


For simplicity, the letters F and g have not been modified when expressed as 
functions of time. For example, the quantities F(t) and g(r) are obtained by replacing 
the angles 6 and <p, respectively, by their appropriate functions of time. The same 
convention applies when the time scale is normalized for F, g, and S. 


A,B,C,D,E coefficients in partial fraction expansion 

a normalized altitude /Radar altitudg \ 

(Radius of planet/ 

b = {l + a 


c speed of light, 2.998 x 10^ meters/second (9.836 x 10^ feet/second) 

F(6>) normalized Muhleman backscatter function (from ref. 2) expressed in 

terms of incidence angle 

F(t) function obtained from F(0) when 8 is related to time 

F(v) function obtained from F(r) by normalizing time scale 

fjj(h) probability density function for h 

fL(i) probability density function for l 

fHL(h,£) joint probability density function for h and l 

G absolute antenna gain 

G 0 maximum antenna gain 

g(cp) normalized antenna gain function 
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g(r) function obtained from g((p) when <fi is related to time 

g(v) function obtained from g(r) by normalizing time scale 

H radar altitude, meters (feet) 

h vertical length of a scattering element 

h(r) received signal resulting from transmitted impulse 

K amplitude factor for pulse echo 

l horizontal length of a scattering element 

M amplitude factor, (1 + a)a;2K 

P(v) pulse response, received signal resulting from transmitted pulse 

P r received power, watts 

Pj- transmitted power, watts 

R radius of planet, meters (feet) 

r distance between radar and arbitrary point on surface, meters (feet) 

S(t) step response, received signal resulting from transmitted step function 

T two-way time delay between radar and closest point on surface, seconds 

t time, seconds 

tp transmitted pulse width, seconds 

u(v) unit step function 

v normalized time variable 

w transformation variable used in integration 



x distance in excess of H between radar and planet along arbitrary ray, 

meters (feet) 

a mean surface slope, <2h/vi 

y planetary angle, degrees 

0 incidence angle, degrees 

X wavelength, centimeters 

p polar coordinate 

a radar cross section, meters^ (feet2) 

0 ^ standard deviation of h 

Oj standard deviation of l 

<Jq radar backscatter function, or radar cross section per unit area 

a Q value of <j q at normal incidence 

r dummy time variable 

Tl time of return from horizon of transmitted impulse 

<p look angle, degrees 

GENERAL DEVELOPMENT 


The general solution for the radar echo step response 


S(v) = K J 


v g2(y) F(v) 


dv 


■X) (1 + v)^ 

is obtained from equation (A21) of the appendix. Equation (A33) 

g 2 (v) dv 


S(v) 


= K f 
J 0 


(i + v)^(i + |)(i + d#] : 
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gives Harrington and Stanley's expression (ref. 1) for the step response in integral form 
(assuming v « 1 and a g 1) where 


and 



H 



The approximate solution to the integral for an omnidirectional antenna (g(v) = 1) is 
given in equation (A34): 


S(v) = 



* + i 1 + 21 + (l+ 


The present development starts with equation (A24) of appendix A: 

1/2 


sin 6 = 


v + 1 


v(v + 2) 


1 + a - 


a 2 v(v + 2) 


from which will be derived an integral for the step response that is less restrictive than 
equation (A33). A solution in closed form is then found for the integral. 

Completion of the multiplication in equation (A24) gives 


sin 6 = — - — /2v + v 2 + 2av + av^ - a^ v^ - a 2 v^ - — — 
1 + vi/ 4 


( 1 ) 


Since cos 8 = /l - sin^e 


cos 8 = — - — ,/l - 2av - av^ + a 2 v 2 + a 2 v^ + 

1 + vV 4 


( 2 ) 


Under the conditions v 2 + 2v = 4 (— ^ a> \ 10" 2 and v 2 + 2v = — , sin 8 can be approxi- 

V a 2 / a 


mated as 


sin e - \/( 2v + v 2 ) (1 ± a) _ b^2v + v 2 

1 + v 1 + v 


( 3 ) 
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and cos 6 can be approximated as 


cos 6 = 


1 

1 + v 


( 4 ) 


where b = ^1 + a. The condition v 2 + 2v = is more restrictive than 

v 2 + 2v S 4 ^ -2 ~ ^ 10 ~ 2 and wil1 be used henceforth. 

Using the foregoing approximations for sin 8 and cos 8, equations (3) and (4) can 
be substituted into equation (A22) for Muhleman’s model to obtain 


F(v) = 


Q! 2 (l + v )2 


(a; + b\/2v + v 2 J* 

where a is the mean surface slope. Substituting equation (5) into equation (A21) gives 
the step response 

S( V ) = Ktt 2 f V ?g 2W dV 

0 (1 + v) ^ + bykv + V 2 ) 


(5) 


( 6 ) 


In equation (6) 


K = 


Go^o 


32tt 2 H 2 (1 + a) 
as in equation (A23) of reference 1. Now, let 

G 0 2 x2ct a-2 

M = (1 + a)a 2 K = 0 0 

32it z H 2 


(7) 


( 8 ) 


and the step response becomes 


S(v) = 


M 

1 + a 


f 

J 0 


Q!g 2 (v) 


(1 + v) (oi + b ^2v + v 2 )‘ 


•dv 


(9) 


As a result, equation (9), the solution in integral form, corresponds to equation (C2) 
of Harrington and Stanley (ref. 1). At this point, their solution is slightly more accurate 
for v < 1 whereas equation (9) is slightly more accurate for v = 1. However, 
Harrington and Stanley continued for g(v) = 1 by finding an approximate solution 
(eq. (C4)) to equation (C2). (See ref. 1.) On the other hand, equation (9) can be solved 
exactly as follows : 
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By letting w = j/2v + v 3 , equation (9) becomes 


S(w) 


m r w 

1 + a *)q 


Q?W 


dw 


|_(1 + w 1 2 )(a + bw) 3 ] 

The integrand can be expanded in partial fractions to give 


( 10 ) 


*r ) . T s-r 

1 + a 


Aw + B 


D 


E 


w 2 + 1 (bw + a) 3 (bw + a) 2 bw + ® 
where the coefficients must satisfy the following set of equations: 


dw 


(ID 


b 3 

0 

0 

0 

b 2 

3ab 2 

b 3 

0 

b 

2ab 

3o! 2 b 

3ab 2 

1 

a 

{a 2 + b2) 

Q! 3 

3a( 2 b 

0 

b 

2ab 

0 

O' 3 

1 

a 

a 2 


( 12 ) 


Equation (11) can now be integrated to give 


S(w) = 


M 

1 + a 


^ log (w2 + 1) + B tan"lw 


C D 

2b(o! + bw) 2 b(a> + bw) 


+ 77 log (bw + a) + — 2_ 
b 2a 2 b 


+ ^ log a 

abb 


(13) 


Solving equation (12) for the coefficients A, B, C, D, and E and substituting equa- 
tion (12) into equation (13) yields the following solution for the integral in equation (10): 


S(w) = 


M 

1 + a 


f 3) log (w2 ♦ 1) ♦ 
2(<* 2 + b 2 ) 3 


ab(3a2 - bg) ^-l 
(o;2 + b 2)3 


w 


a c 


oja? - b2) 


2(q; 2 + b 2 )(a + bw)^ (a 2 + b 2 ) 2 (a + bw) 


a?2(q;2 _ g^2) 

(a2 + b2) 3 


log (a + bw) 


1 Qf 2 - b2 + q? 2 (q? 2 _ 3b 2) ^ 

2(a; 2 + b 2 ) (a 2 + b 2) 2 (a 2 + b 2 ) 3 


(14) 
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The solution for the step response in terms of normalized time and altitude can now be 
found by rearranging equation (14) and transforming back to the original variables to give 
the final result: 


S(v) = 


M 

1 + a 


c 

q'2(q; 2 - 3a - 3) , J a; 2 (v 2 + 2v + 1) 

o -xS log \r , =1 


2(cx 2 + a + l)^ 


a + l/(2v + v 2 )(l + a) 


+ a Vl + a (3o? 2 - a - 1) tan _ i 

(cx^ + a + 1)^ 



cr 


2 (cx 2 + a + 1) [cx + ^(2v + v2)(l + a)]' 

cx(cx 2 - a - 1) 

(cx2 + a + l)^ |cx + ^/(2v + v 2 )(l + a)] 2(cx2 + a + 1)' 


3a 2 - a - 1 


\ 


(15) 


As in reference 1, the pulse response can now be found by taking the difference of two step 
responses, that is, 


P(v) = S(v) u(v) - S (v - -£) u (v - ^ 


(16) 


where tp is the transmitted pulse width. 


Restrictions 

As previously stated, Harrington and Stanley (ref. 1) assumed a = 1 and v « 1 

i n- 2 

(e.g., v = 10-2). The solution in equation (15) assumes v 2 + 2v = — — . These restric- 
tions are plotted in figure 1 in terms of the normalized variables a and v to show the 
domains in which the two solutions are valid. In figure 2 the restrictions are plotted in 
terms of the variables t (time) and H (altitude) for a specific case, namely the planet 
Mars. As shown by figure 2, the solution of equation (15) is valid for long pulse widths 
whereas the previous solution of Harrington and Stanley (ref. 1) is not. 


RESULTS 


In order to demonstrate the utility of the present solution and to provide examples 
for comparison between this solution and the previous one of Harrington and Stanley 
(ref. 1), the expected values of the received pulse amplitude as a function of time were 
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calculated for several values of altitude H and the mean surface slope a. The results 
are plotted in figures 3 to 5. (A rectangular transmitted pulse of 10-psec duration and 
a planet radius of 3370 km (~ 1.1 X 10^ ft) were used in all cases, and square-law detec- 
tion was assumed.) The figures show the results obtained using the present solution 
(eq. (15)), using the previous solution of Harrington and Stanley (eq. (A34)), and using 
numerical integration of the exact integral solution (eq. (A29)). (Integration intervals 
of 2.5 nsec for a = 1.0 and 1.25 nsec for a = 0.01 were used in the numerical integra- 
tion.) Note that the planet radius used is that of Mars; thus, the restrictions on time and 
altitude shown in figure 2 are applicable to the examples in figures 3 to 5. 

The results for an altitude of 1.52 km (5000 ft) and for a of 1.0 and 0.01 are 
shown in figures 3(a) and 3(b). Good agreement exists between the results of the present 
solution and the results of numerical integration, the only discernible difference occurring 
near the pulse peak for a = 0.01 (fig. 3(b)). However, the previous results of Harrington 
and Stanley (ref. 1) are quite different in both cases. It can be seen from figure 2 that the 
previous restriction on time for H = 1.52 km (5000 ft) is exceeded at t = 0.1 psec, or 
long before the pulse peak occurs. It can be seen from figure 3(b) that the previous solu- 
tion begins to decrease for time less than one pulse width, but such a decrease does not 
occur physically. This anomaly accounts for the trailing-edge amplitude becoming nega- 
tive in that case. 

The results in figure 3(c) are for an altitude of 152 km (500 000 ft) and for a = 0.01. 
In this case, Harrington and Stanley's previous restriction on time is not exceeded until 
t = 10 psec; and their result closely agrees with the results of the other two solutions. 

The amplitudes in figure 3 have all been normalized by dividing by the factor M. 
Therefore, comparisons of peak amplitude, as well as pulse shape, are valid among the 
three solutions for a given a and H. 

In figure 4 the results for a - 1.0 and H = 1.52 km (5000 ft) are plotted after the 
amplitude at the end of one pulse duration has been normalized to unity. Such a presenta- 
tion emphasizes the differences in pulse shape among the three solutions. 

The solutions for the same return pulse are plotted in figure 5, in this case after 
passage through a third-order low-pass Butterworth filter with a 3-dB cutoff frequency 
equal to 0.5/tp. The filtering was accomplished using the technique described by 
Harrington and Stanley in reference 1. A slight difference between the previous solution 
of Harrington and Stanley and the other two solutions can again be seen. 

The advantage of the present solution over the previous solution of Harrington and 
Stanley is improved accuracy for cases where their restrictions on time and amplitude are 
violated. The advantage of the present solution over the solution by numerical integration 
is a reduction in computer time. By way of illustration, the computer central processing 
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unit time required at H = 1.52 km (5000 ft) and a = 0.01 was 4.9 seconds for the solu- 
tion by numerical integration compared to 0.15 second for the present solution. 

CONCLUSIONS 

The solution presented can be used to compute the expected return pulse amplitude 
and shape for planetary radar altimeters, where the surface backscatter characteristics 
are described by Muhleman’s model. The solution is valid over a wide range of altitude 
and pulse widths, and results obtained with this solution show excellent agreement with 
results obtained by numerically integrating the exact integral solution. Computations 
using this solution require significantly less computer time than computations using 
numerical integration. In addition, the solution is less restrictive on the normalized alti- 
tude and time than the solution previously presented in NASA TN D-5220 by Harrington and 
Stanley. Furthermore, the new solution is more accurate at low altitudes and for long 
pulse widths. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 28, 1971. 
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APPENDIX A 


DERIVATION OF RADAR ECHO STEP RESPONSE 


The derivation of the radar echo step response in this appendix is a slight modifica- 
tion of the one presented by Harrington and Stanley in reference 1. Since the electro- 
magnetic wave transmission properties of a vacuum and of the atmosphere and the reflec- 
tion characteristics of a planetary surface are independent of power (for reasonable power 
levels), the atmosphere and surface can be treated as a linear system, and linear system 
theory can be applied to the problem of finding the power received at the altimeter 
(system output) in response to a transmitted pulse (system input). The first step is to 
find the system impulse response, that is, the received power as a function of time due 
to an impulse transmitted at time t = 0. Note that in this case the input and output 
functions are in units of power rather than amplitude, as is the more usual case. 

From the fundamental radar range equation the received power is given by 


P t G2\2 g 

r (47r) 3 r 4 


(Al) 


For the radar altimeter the geometry of the model is shown in figure 6. Consider a unit 
impulse transmitted at time t = 0. The energy reflected by the annular ring dA, which 
is at a constant range r = H + x from the altimeter, will be received by the altimeter at 
time t = T + t, where 


T _ 2H 
c 

(A2) 

11 

(A3) 


The magnitude of the energy received from dA is 


h(r) dr = 


G 2 A 2 cr 0 dA 
(4ir) 3 (H + x) 4 


G 2 A 2 a 0 dA 

<4ir)SH4(l + i) 4 


Now the antenna gain is a function of the look angle <p and can be expressed as 


(A4) 


G = G og (<#>) 


(A5) 


where 


0 § g(<f>) § 1 


11 



APPENDIX A - Continued 


For an isotropic surface the backscatter function a 0 is a function of the incidence 
angle 6, but not of the azimuth angle. Thus the function a Q can be expressed as 

CT o = ^o F ( 0 ) (A6) 

where a Q is the backscatter per unit area at normal incidence. The differential 
area dA is 


dA = 2rR 2 sin y d y 


(A7) 


Combining equations (A4) to (A7) results in 


h(r) dr = 


Go 2 g 2 (<£) r 2 s i n y d y 


32ir 2 H 4 fl + £ 


4 


Applying the law of cosines to triangle POQ in figure 6 gives 

x 2 + 2xH 

2R(H + R) 

By differentiating equation (A9), 

H 2(i + i.\ dT 

sinydy = ( x ^ H) dx = V TJ 
R(H + R) TR(H + R) 


(A 8) 


(A9) 


(A10) 


is obtained. Substituting equation (A10) into equation (A8), and letting 


a = H 

R 


(All) 


and 


K 


Gq 2 ^- 2 Uq 

32tt 2 H 2 (1 + a) 


(A12) 


gives finally 


h(r) dr = 


Kg 2 (c/>) F (0) dr 


Til + 


if 


(A13) 


Dividing by dT and expressing g(<p) and F(0) as functions of r gives the following 
system impulse response which corresponds to equation (A24) of reference 1: 
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APPENDIX A - Continued 


Kg 2 (T) F(t) 


h(r) = \ 


^ + t)' 


(0 S T S t l ) 
(t < 0, t > r L ^ 


(A14) 


where 


t l = T -1 + Jl+f 


(A15) 


The limit r = 0 is the time when the transmitted impulse first returns from the surface. 
The limit corresponds to reflection at the horizon (0 = 90°). 

From linear-system theory, the power received from an arbitrary transmitted sig- 
nal Pt(t) can be found by convolution as 


P r (t) = J h(r)Pt(t - r)dT 

_ oo 

Upon substitution of equation (A14) this becomes 

P r (t) = f - r) dr 

° ('^) 

To find the step response S(t), let Pt(t) = u(t). Then 


(A16) 


(A17) 


S(t) = f f ^ g2(T) F(T) u(t - T)dT 
T J 0 / t\3 


(A18) 


>4 


Equation (A18) can be written as follows: 


S(t) = / 


r 

K 

§ 2 (t)f(t) dT 

T ^ 

H)‘ 

K I 

r T L g 2 (r)F(r) 

T . 



(‘ s TL ) 


(t > TL ) 


(A19) 
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APPENDIX A - Continued 


Now, introduce the normalized time variable v, where 


Then 



x_ 

H 


r K r T g2( v) F ( v) dv 

J 0 (1 + v) 3 


S(v) = < 


K 


f T l/ T g 2 (v) F(v) dv 

(1 + v) 3 



Thus far only backscatter functions which are an arbitrary function of the 
angle 8 have been considered. Henceforth, attention will be restricted to the 
model (ref. 2) for F (8), namely 


F(0) = 


a 3 cos 8 

o 

(sin 8 + a cos 8)° 


Applying the law of sines to triangle QOP in figure 6 gives the relationship 


sin y _ sin 6 
H + x h + R 


or 


sin 8 = / + a ■ sin y 
a(l + v) 


Using the identity sin y = fL - cos 2 r and equation (A9) results in 


sin 8 = 


1 + v 


r(v + 2) 


o -Tl !/ 2 


which corresponds to equation (B5) of reference 1. 
Now let 

r i— 

>2. 


A(v) = <v(v + 2) 


-iV / 2 


1 + a - ]j> 


(A20) 


(A21) 


incidence 

Muhleman 


(A22) 


(A23) 


(A24) 


(A25) 
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APPENDIX A — Continued 



so that 


sin 


e = A (zl 

1 + V 


(A26) 


and 


cos 8 = 


T i-|((l + v)2-A2(v)»^i 

J + V ' 1 + V 


(A27) 


Then Muhleman's model becomes 


F( v ) = o' 3 (l + v) 2 B(v) 

|a(v) + <*B(v)] 

By substituting equation (A28) for F(v) in the integral of equation (A21), the step 
response becomes 

o/ \ ^ 3 f V g 2 (v) B(v) 

® (1 + v)[a(v) + ci!B(v)] 


(A28) 


(A29) 


Thus far no approximations have been made. However, the integral in equation (A29) 
is quite unwieldy and has not been integrated in closed form. It could, of course, be inte- 
grated using numerical-integration techniques, and this was done to check the validity of 
the approximate solutions (as discussed in the main body of this report). 

At this point, Harrington and Stanley (ref. 1) made some assumptions in order to 
obtain an analytic approximation to the step response. With the assumptions that a = 1 
and v « 1, Harrington and Stanley approximated sin 6 by equation (B7) of reference 1 
as 


sin 9 



1 + v 


(A30) 


and cos 6 by equation (B6) of reference 1 as 


cos 6 ~ 


1 

i + v 


(A31) 


With these approximations the Muhleman backscatter function becomes (eq. (B8) of ref. 1) 



(A32) 
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APPENDIX A - Concluded 


Substitution of equation (A32) into equation (A21) gives equation (C2) of reference 1 
for the step response as follows: 



g 2 (v ) dv 


( 1+V )^( 1+ !( 1+ ^f 


(A33) 


For an omnidirectional antenna (g(v) = 1), the step response integral (eq. (A33)) was 
approximated in equation (C4) of reference 1 by the following analytic expression: 


S(v) ~ 


Ka; 2 Jl 1 

_ I _ 

1 


T 

2 

[♦MW 


(A 34) 
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APPENDIX B 


REVIEW OF MUHLEMAN'S MODEL FOR THE 
RADAR BACKSCATTER FUNCTION 


In reference 2, Muhleman derived an expression for the radar backscatter function 
from a random surface using the principles of geometric optics. Muhleman assumed that 
the surface consisted of individual plane scattering elements, all of which are long com- 
pared with the wavelength of the incident radiation. By geometric optics a scatterer will 
reflect in a given direction only if the normal to the scattering surface is coplanar with the 
incident and reflected rays and if the normal bisects the angle between these rays. For 
radar backscatter this amounts to the scattering element normal being parallel to the inci- 
dent ray. Muhleman utilized the assumed statistical properties of the surface scattering 
elements to determine the probability of the normal meeting this requirement. In his 
derivation Muhleman assumed an isotropic surface, and he did not consider the effects of 
polarization. 

The surface is characterized by length l parallel to the mean surface and a height 
h perpendicular to the mean surface. If h and l are assumed to be independent ran- 
dom variables, then joint probability density function fHL(h,£) can be expressed as the 
product of the individual density functions fH(h) and fL,(Z): 


fHLM = f H (h) f L G) 


(Bl) 


By transforming to polar coordinates p and 6 and integrating over p, the probability 
density function for the normal can be obtained in terms of 6, the angle between the scat- 
tering element normal and the mean surface normal. The angle 6 is also the angle of 
incidence relative to the mean surface. 

Muhleman (ref. 2) derived the backscatter function for two cases. In one case he 
assumed a Gaussian distribution for the height variable h and a Rayleigh distribution for 
the length variable l. In the other case, he used exponential and Poisson distributions for 
the height and length variables, respectively. The resulting scattering functions for both 
cases were fitted by Muhleman to existing radar backscatter data from the Moon and 
Venus. Although both models exhibited good agreement with the data, the exponential 
model produced the best fit; hence, it was utilized by Harrington and Stanley (ref. 1) and in 
the work described in this report. 

Muhleman’s model for the backscatter function derived for the exponential case is 


F (6) = 


af3 cos 0 

(sin 9+oi cos 6)^ 


(B2) 
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APPENDIX B - Concluded 


where 



a 

mean surface slope, 

°h/°i 

a h 

standard deviation of 

h 


standard deviation of 

l 


The function has been normalized to unity at vertical incidence (8 = 0); thus, the radar 
cross section per unit area can now be expressed as 

a o (e) = 6 o F(0) (B3) 

where d 0 is the cross section per unit area at vertical incidence. 

In comparing his theoretical model to the radar reflection data for the Moon, 
Muhleman (ref. 2) adjusted the value of the parameter a to obtain the best fit between 
his model and the measured data. At wavelengths of 68 cm, 10 cm, and 3.6 cm, he 
obtained values for a of 0.145, 0.28, and 0.35, respectively. Using data from radar 
observations of Venus, Muhleman estimated that the value of a lies between 0.04 and 0.12 
for that planet at a wavelength of 68 cm. In order to span the values of a obtained by 
Muhleman (ref. 2) and to facilitate comparison between Harrington and Stanley’s results 
(ref. 1) and the results contained in this report, values for a of 0.01 and 1.0 were 
selected. 
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Figure 1.- Restrictions on normalized time and altitude. 
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(a) a = 1.0; H = 1.52 km (5000 ft). 

Figure 3.- Comparison of return pulse shapes and amplitudes. 
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(b) a = 0.01; H = 1.52 km (5000 ft) 
Figure 3.- Continued. 
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Figure 5.- Comparison of filtered pulses 
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